Fractal dimension of domain walls in two-dimensional Ising spin glasses 
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We study domain walls in 2d Ising spin glasses in terms of a minimum-weight path problem. 
Using this approach, large systems can be treated exactly. Our focus is on the fractal dimension 
df of domain walls, which describes via {£) ~ L t the growth of the average domain- wall length 
with systems size L x L. Exploring systems up to L = 320 we yield df = 1.274(2) for the case 
of Gaussian disorder, i.e. a much higher accuracy compared to previous studies. For the case of 
bimodal disorder, where many equivalent domain walls exist due to the degeneracy of this model, 
we obtain a true lower bound df = 1.095(2) and a (lower) estimate df = 1.395(3) as upper bound. 
Furthermore, we study the distributions of the domain-wall lengths. Their scaling with system size 
can be described also only by the exponent df, i.e. the distributions are monofractal. Finally, we 
investigate the growth of the domain- wall width with system size ("roughness") and find a linear 
behavior. 
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I. INTRODUCTION 

In this paper we examine the scaling behavior of 
minimum-energy domain wall (MEDW) excitations for 
two-dimensional Ising spin glasses. Spin glasses are a 
prototypical model [JH, d, BJ for systems with quenched 
disorder in statistical mechanics. In general, despite more 
than two decades of intensive research, many properties 
of spin glasses, especially in finite dimensions, are still 
not well understood. The situation is slightly better for 
two-dimensional spin glasses, because it is now widely 
accepted that no ordered phase for finite temperatures 
exists d, 0, 0, HI, 9- For d = 2 the behavior can be de- 
scribed well by a zero-temperature droplet scaling (DS) 
approach [ljj [ll|, QjJ ■ Within this approach, the excita- 
tion energy of domain walls and other excitations scales 
like AE ~ L e , where L is the system size and 9 is re- 
ferred to as stiffness exponent. Within DS, 9 is assumed 
to be universal for all types of excitations. Further, the 
surface of an excitation is assumed to display a scaling 
that is characterized by a fractal dimension df. 

For Gaussian disorder of the interactions, prior inves- 
tigations of domain walls resulted in estimates for the 
stiffness exponent 9 = —0.287(4) 0, [H| and the values 
for the MEDW fractal dimension listed in Tab. E For 
this model, it was recently indeed confirmed [3, Ell that 
the value of 9 is the same also for other types of excita- 
tions. Recent studies suggested that such domain walls 
are governed by stochastic Loewner evolution (SLE) pro- 
cesses [ill [TtJ ■ Within conformal- field theory, it seems 
possible to relate the MEDW fractal dimension df to the 
stiffness exponent 9 via 



d/-l = 3/[4(3 + 0)]. 



(1) 



high accuracy is needed to verify whether the proposed 
relation is compatible with the data. Note that the er- 
ror bars of the previous results for df are typically of 
order 10 -2 or larger, i.e. ten times larger than our high- 
precision result. Also, we reach much larger system sizes 
than previous work (also a bit larger than the recent 
study in Ref. [l8|], which is anyway for a different lat- 
tice type), which also reduces systematic errors due to 
unknown corrections to scaling. 



Reference 


df 


Geom. 


System 


Alg. 


Middleton [19] 


1.25(1) 


tr 


256 x 256 


M 


Bray/Moore [20] 


1.26(3) 


sq 


12 x 13 


TM 


Kawashima/Aoki [21] 


1.28(2) 


sq 


48 x 48 


M 


Bernard et. al. [171 


1.28(1) 


tr 


720 x 360 


M 


Rieger et. al. [5J 


1.34(10) 


sq 


30 x 30 


BC 


Palassini/ Young [221 


1.30(8) 


sq 


30 x 30 


BC 


Weigel/ Johnston [18] 


1.273(3) 


hex 


256 x 256 


M 



Note that for a wide range of values of 9, one gets very 
similar results for the fractal dimension, e.g. df(9 = 
-0.2) = 1.268, while d f (6 = -0.3) = 1.278. Hence, a 



TABLE I: Previous results on the fractal dimension of 
MEDWs arising from Gaussian disorder in two dimensions. 
From left to right: Reference, estimate value of the fractal 
dimension df, geometry of the system (sq: square lattice, 
tr: triangular lattice, hex: hexagonal lattice), largest system 
studied and numerical algorithm used, for details see cited ref- 
erences (M = matching approaches, TM = (T = 0) transfer 
matrix method, BC = branch-and-cut algorithm). 



For the ± J-model it was found that the MEDW en- 
er gy s aturates at a nonzero value for large system sizes 
[JUil, which means 9 = 0. It exhibits a high degener- 
acy of ground states (GSs) and thus allows for numer- 
ous MEDWs, e.g. varying in length, see Fig. U As a 
result, the concept of a MEDW is not clear-cut. Refer- 
ring to this model, the SLE scaling relation above can- 
not be adopted. Recent attempts to capture the frac- 
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FIG. 1: For the ±J model, a high degeneracy appears. 
Hence, for each of the GSs, many MEDWs separating the 
left and right border exist, which having all the same min- 
imum energy. Left: many MEDWS for a sample system of 
size L = 80, middle (right): corresponding ±J-min (±J- 
max) MEDW. Dashed lines denote free boundary conditions, 
solid lines indicate periodic boundary conditions. 



tal properties of typical MEDWs arising from bimodal 
disorder concluded with values df — 1.30(1) [H| and 
df = 1.283(11) [Hj]. Nevertheless, in these works there is 
still no systematics concerning the selection of a represen- 
tative MEDW, as mentioned in Ref. [l||. The sampling 
of the domain walls is not controlled, hence different con- 
figurations having the same energy do not contribute to 
the results with the same weight. Currently, no fast al- 
gorithm is known which allows to sample the degenerate 
GSs and/or MEDWs correctly. For a correct sampling, 
all degenerate configurations must contribute to the re- 
sults with the same weight or probability. So far this can 
be done only for small systems through enumeration of all 
GS configurations [25| . Some investigations use a scaling 
relation from droplet theory [ll[, which quotes that the 
variance of the DW entropy should also scale like a power 
law with the system size, with the same exponent df. In a 
recent investigation based on this relation for square lat- 
tices of size up to L = 256 [26[ a value of df = 1.090(8) 
for 24 < L < 96 and d f = 1.30(3) for L > 96 was re- 
ported. On the other hand, other studies on the scaling 
behavior of the MEDW entropy [13, HH suggested that 
the behavior of zero-energy DWs is not consistent with 
the droplet scaling picture. In this context it was pro- 
posed to treat MEDWs of zero and non zero energies as 
distinct classes. These findings were supported by calcu- 
lations for square lattices as well as for an aspect ratio 
different from unity. Hence, due to these results, a direct 
determination of the MEDW fractal dimension appears 
to be preferable. 

So as to shed further light on their fractal properties, 
we investigate MEDWs originating from bimodal and 
Gaussian disorder. We use exact algorithms [2{|, which 
allow to obtain the domain walls with the true lowest en- 
ergy. Here, we describe a novel approach to the problem 
of finding MEDWs in terms of minimum- weight paths, 
detailed in the forthcoming section. It allows to put a 
lower bound and an upper estimate on the scaling behav- 
ior of MEDWs and it might further provide a sound basis 
to investigate the scaling behavior of typical MEDWs for 
the ±J model, similar to the random-field Ising model, 
where all ground states can be represented by a certain 



graph [30(, which should allow for unbiased sampling of 
ground states. Moreover, the corresponding picture of 
minimum-weight paths on undirected graphs, that allow 
for negative edge weights, leads to a percolation problem 
that appears to be interesting on its own, i.e. seems to be 
in a new universality class [3 11 ] . For the bimodal system, 
to avoid aforementioned problems of sampling correctly 
equivalent MEDWs of different length, we ask for those 
MEDWs with extremal lengths, bearing a true lower 
bound and an upper estimate 1.095(2) < df < 1.395(3) 
in case of bimodal disorder. In addition, for the Gaus- 
sian system, we study square lattices with a large num- 
ber of samples and much larger system sizes than pre- 
viously studied in the literature, yielding a result for 
the fractal dimension df = 1.274(2), which has an en- 
hanced precision compared to the results shown in Tab. 
Uand compares well with the value e^ LE = 1.2764(4) ob- 
tained via inserting 9 = -0.287(4) into Eq. (T]). Finally, 
we also study for the first time the scaling-behavior of 
the domain-wall roughness and of the distributions of 
domain- wall length and width. 

The paper is organized as follows. In the next section, 
we present the details of the model and the numerical 
algorithms we have applied. In the third section, we show 
all our results. We conclude with a summary in section 
four. 



II. MODEL AND METHOD 

In the framework of this paper, we present GS calcula- 
tions of two-dimensional Ising spin glasses with nearest- 
neighbor interactions. The model consists of TV = L x L 
spins a = (01, ... , ctn) with cr,; = ±1 located on the sites 
of a square lattice. The energy is given by the Edwards- 
Anderson Hamiltonian 

H {a) = - Jij (2) 

where the sum runs over all pairs of neighboring 
spins, with periodic boundary conditions (BCs) in the 
x— direction and open BCs in the y— direction. The 
bonds Jij arc quenched random variables drawn accord- 
ing to a given disorder distribution. They can take either 
sign and thus lead to competing interactions among the 
spins. Here, we consider two kinds of disorder distribu- 
tions: (i) Gaussian with zero mean and variance one, and 
(ii) a bimodal distribution = ±1 with equal probabil- 
ity (±J model). 

In this work, we study minimum-energy domain walls, 
which are certain excitations that are defined, for each 
realization of the disorder, relative to two spin config- 
urations: cr is a ground state with respect to periodic 
BCs (a;— direction) and cr AP , a GS obtained for antiperi- 
odic BCs. One can realize antiperiodic BCs by invert- 
ing the signs of one column of spins, described by the 
Hamiltonian H AP (a). Considering both spin configura- 
tions, MEDWs separate two regions on the spin lattice: 
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one, where the spins have the same orientation in both 
GSs and another, where the spins have different orien- 
tations regarding the two GSs. fn the above periodic- 
antiperiodic setup, one domain wall will run through the 
inverted bonds. This straight-line domain wall is not of 
interest to us. Instead, we are interested in the MEDW 
which consists of those bonds that arc fulfilled in exactly 
one of the two GSs. Hence, the energy of the MEDW 
is given by H AP (a AP ) — H(a) and this is the minimum 
energy among all the system-spanning (top-down) DWs. 
This MEDW feature is an integral part of our attempt 
to study the problem of finding MEDWs in terms of a 
minimum- weight path problem. 

First, we will now outline the main steps of our al- 
gorithm and elaborate on them afterwards. The algo- 
rithm can be decomposed into the following 4 steps: For 
each given realization of the bond disorder (i) find a 
GS spin configuration consistent with periodic BCs in 
x— direction, (ii) determine all possible MEDW segments 
on the dual of the spin lattice, (iii) map the problem to 
an auxiliary graph and find a minimum weighted per- 
fect matching (MWPM), (iv) interpret the MWPM as 
minimum-weight path on the dual graph that represents 
a MEDW on the spin lattice. 

Step (i): For the calculation of the GS a. we apply 
an exact matching algorithm that works for planar sys- 
tems without external fields, e.g. a square lattice with 
periodic BCs in at most one direction. For this purpose, 
the system has to be represented by its frustrated pla- 
quettes and paths connecting those pairwise. Finding 
a minimum-weighted perfect matching on the graph of 
frustrated plaquettes then corresponds to finding a GS 
spin configuration on the spin lattice. For a compre- 
hensive description, see Refs. [29|, |H, [H, This well 
established method yields the exact GS of the frustrated 
spin-glass model and allows to explore large system sizes, 
easily more than 10 5 spins, within a reasonable amount 
of computing time. 

Step (ii): Once the GS is obtained we construct the 
dual of the spin lattice as weighted graph G = (V,E,u>). 
The vertices V are given by all the distinct plaquettes on 
the spin lattice and edges e£E connect vertices, where 
the corresponding plaquettes have a bond in common, i.e. 
the edge crosses the corresponding bond, see Fig. [21 Note 
that there are two "extra" vertices above and below the 
system. Further, a weight (or distance) u>(e) is assigned 
to each edge of G, equal to the amount of energy that 
it would contribute to a MEDW, i.e. u>(e) = 2Jij<Ji<jj 
for being the bond crossed by the edge e, a being 
the GS obtained in step (i). In this sense, the dual graph 
comprises all possible MEDW segments. A DW is a path 
in the dual graph and the energy of a DW is the sum of 
the weight of all segments being part of a DW. The GS 
property of the spin configuration has an impact on the 
weight distribution of G. There are negative edge weights 
but one cannot identify loops with negative weight. This 
will be of importance later. In summary, a MEDW is a 
path with minimum weight joining the extra vertices of 
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FIG. 2: (color online) Construction of the dual of the spin 
lattice. Left: Spin lattice with periodic BCs in the horizontal 
direction (denoted by black bonds at the left and right col- 
umn), and free BCs along the top and bottom boundary. The 
sites are denoted by the circles, the bonds are the straight and 
jagged lines. For clarity, spins are not drawn. Right: Vertices 
and edges of the dual graph. Note that two extra vertices on 
top and at the bottom are introduced to account for the free 
BCs. 



G. 

Step (iii): Since G is an undirected graph that allows 
for negative edge weights, it requires matching techniques 
to construct minimum- weight paths 35[. Following this 
reference, we therefore map the problem onto an auxiliary 
graph G a obtained from the dual graph by the following 
procedure: For every vertex i in the dual graph (except 
the two extra vertices above and below the lattice) Ga 
contains a pair of "duplicate" vertices i^ a \ which are 
connected by an edge of zero weight. Furthermore for 
every edge e = (i,j) in the dual graph, two additional 
vertices ay'" , b^'^ are introduced which are connected 
by an edge (ay'^, b™>). One of these two vertices is 
connected to the duplicate vertices representing i (e.g. 
via edges (a^^jW), (a^'i\i^)), the other one to the 
duplicate vertices representing j (via edges (b^'^ 

j( b >)). All edge weights are zero, except for the 
edges connecting one of the additional vertices to the 
duplicates it is connected to, which carry the weight of 
uj(e), see Fig. EH We then apply an algorithm from the 
LEDA HIES} libr ary to obtain a MWPM on Ga- 

Step (iv): The MWPM consists of a certain subset of 
the edges of Ga- In order to interpret it as minimum- 
weight path on the dual graph, one has [35[ to perform a 
partial inverse transformation of step (iii). This means, 
we merge all adjacent vertices with the same type, i.e. 
all pairs i^ a >,i^ of "duplicate" vertices and all pairs 
a 0>j) ( 5OJ) f additional vertices. The edges in the match- 
ing which are between vertices of the same type disap- 
pear, while the other edges "remain" , which means each 
edge (a/b( t 'i\i/j( a W( b ') becomes an edge (a/b^\i/j). 
Thus, some of the edges contained in the matching dis- 
appear, while the remai ning ones form (for a proof and 
further details see Ref. [351 ) a path of minimal weight, 
connecting the extra vertices of G. This path in turn 
corresponds to the MEDW on the spin lattice. 

One advantage of the above procedure, in comparison 
with computing the GSs for periodic and antipcriodic 
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FIG. 3: Construction of the auxiliary graph: Every connected 
pair of vertices on the dual is replaced by a construct of 6 
vertices and 7 edges. This transformation is illustrated for 4 
edges depicted on the left, with weights wl, w2, w3 and w4. 
After the transformation the gray edges carry zero weight and 
the black edges carry the same weight as the original edge on 
the dual. For an illustrative purpose the vertices are divided 
into round and squared ones. 




FIG. 4: Extension of the algorithm: subpath ci, replacement 
path C2 and pivoting spin cluster S that might be flipped in 
order to improve a MEDW regarding its length. 



boundary conditions, is that it yields an explicit repre- 
sentation of a MEDW. This allows to determine quanti- 
ties like length/roughness scaling exponents directly. The 
main advantage is related to the ± J-model. There, it is 
possible to lift the degeneracy among the MEDWs. One 
can obtain MEDWs with an exactly minimum and a max- 
imal (i.e. not maximum) number of segments, i.e. exactly 
shortest and very long MEDWs. 

First, the slight modification u)(e) — > ui(e) + e (e > 0) 
represents a negative feedback for the inclusion of path 
segments, hence the MEDW will be among all MEDWs 
that one which includes a minimal number of path seg- 
ments, i.e. a shortest MEDW. The value of e must 
be small enough to maintain the minimum-weight path 
structure on the dual graph, e.g. e < 1/\E\. This mod- 
ification is referred to as ±J-min. The modification 
u)(e) — > to(e) — e represents a positive feedback for the 
inclusion of path segments to the MEDW. Hence, the 
MEDW obtained in this way will be, among all MEDWs, 
one which will include a large number of path segments. 
This is referred to as ±J-max. However, in the latter 
case arises a serious difficulty: the weight distribution 
may allow for loops with negative weights now, and, due 
to the nature of the MWPM problem, the algorithm re- 
turns a minimum-weight path in the presence of loops 
with negative weight. Hence, only the total number of 



segments of the minimum-weight path together with all 
loops (of zero energy in the unmodified graph) is maxi- 
mized, not the number of segments of a minimum- weight 
path alone. So we are only able to obtain a lower bound 
on the longest MEDWs via using the ± J-max approach. 
Note that obtaining the true longest minimum-weight 
MEDWs is an NP-hard problem, which means that only 
algorithms arc known, where the running time increases 
exponentially with the number N of spins. 

Another drawback is the fact, that the MEDW de- 
pends on the spin configuration determined in step (i). 
To get really the shortest MEDWs (± J-min2), we there- 
fore allow the algorithm to change the GS, if this leads to 
a shortening of the MEDW. The basic idea is that differ- 
ent GSs differ by a finite number of zero-energy clusters of 
reversed spins. We are only interested in clusters, where 
part of the boundary coincides with a MEDW, which 
lead to a different MEDW of the same energy when flip- 
ping the cluster. Note that also the edges connecting to 
the two extra nodes are considered here as part of the 
MEDW, hence a flip of one or several zero-energy clus- 
ters might lead to another MEDW which has nothing in 
common with the starting MEDW. Technically, we look 
for the shortest MEDWs by finding replacement paths 
for certain subpaths of the MEDW, depicted in Fig. |U 
Firstly, one has to find a subpath ci of the MEDW with 
weight lu{c\) < 0. If there is an replacement path C2 with 
uj{c2) = |w(ci)| and if len(c2) < len(ci) flip the cluster of 
spins surrounded by c\ and ci to yield a shorter path with 
minimum weight. This is repeated until no such shorten- 
ing of the MEDW is possible any more. Since the flipping 
of the spin cluster does not cost energy (the loops ci , C2 
are zero-energy loop in the unmodified graph), the GS 
property of the spin configuration is maintained. This 
latter elaboration of the algorithm is computationally 
more expensive because we have to consider all possible 
subpaths ci with uj(c\) < 0. This means, we can obtain 
the true shortest MEDWs, considering all possible GSs, 
only for small sizes L < 32. For large sizes, we study 
only the shortest MEDW for the given GS, which was 
obtained in step (i). 

III. RESULTS 

In Fig. [1] we show several MEDWs for a system with 
±J disorder, one long ±J-max MEDW and one short- 
est ±J-min MEDW, respectively. Already from these 
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TABLE II: Number of samples investigated for the different 
system sizes. 
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FIG. 5: (color online) Scale dependence of the average MEDW 
length. Fits to the form {£) ~ L d i (dashed lines) yield the 
exponents quoted in Tab. lIIII The inset shows the local slopes 
for the Gaussian and ± J-min(max) MEDWs as function of 
1/L, providing comparable values for df but with larger error 
bars (see text). 



samples one can expect an observable difference in their 
scaling behavior. 

We expect the disorder-averaged MEDW length £ to 
scale with the system size according to (£) ~ L df , where 
(. . .) denotes the disorder average and 1 < df < 2 de- 
notes the fractal dimension of the MEDW. Its roughness 
r, given by the difference in the coordinates of its left- 
most and rightmost position on the spin lattice, should 
display the scaling behavior (r) ~ L dr , with a rough- 
ness exponent d r . To determine the scaling behavior of 
the MEDWs, we have studied systems with sizes up to 
L = 320 and averaged over up to 40000 realizations of the 
disorder, see Tab. HU 

In Fig. [5l we show the result for (£) as function of 
system size for the Gaussian disorder and the three cases 
studied for the ± J model. The data can be fitted very 
well to power laws, the results are shown in Tab. IIIIl As 
an alternative method, we also estimate the value of the 
scaling exponent df using the local (successive) slopes of 
the data points, see also [39|]. The results are compatible 
with the data given in Tab. lIIll only this procedure leads 
to a more conservative estimate of the numerical error, 
leading to the error bars as given in Tab. IIIIl 





df 


Q 


d r 


Q 


Gaussian 


1.274(2) 


0.45 


1.008(11) 


0.40 


±.I mill 


1.095(2) 


0.27 


1.006(6) 


0.78 


±J-max 


1.395(3) 


0.16 


0.993(8) 


0.35 



TABLE III: Scaling exponents of the mean MEDW length 
(df) and its roughness (d r ). Fits are restricted to L > 26 
(df) and L > 50 (dr). The value of Q gives a measure for the 
quality of the fit [3g| . 



The estimate for the MEDW fractal dimension in case 
of a Gaussian distribution of the interaction strengths is 
in good agreement with earlier results, but has an en- 
hanced precision, see Tab. HI Note that the treatment 
of large systems reduces the influence of systematic er- 
rors due to unknown corrections to scaling, hence pro- 
viding a very reliable result for df. Further, it is con- 
sistent with the result dj LE = 1.2764(4) predicted by the 
SLE scaling relation, where we considered 8 = — 0.287(4) 
from Ref. [13]. For the ±J-model we obtained a lower 
bound df = 1.095(2) that is distinct from 1, indicating 
that overhangs are still significant for MEDWs with min- 
imal length. Further, the estimate of the upper bound 
df = 1.395(3) points out that DWs with maximal length 
are not space-filling. Note that using the past results 
[3 HI| , which are based on an uncontrolled sampling of 
domain walls, one could not exclude these values df = 1 
and df = 2. 

As pointed out in the description of the algorithm, the 
± J-min value actually overestimates the scaling behav- 
ior of MEDWs with minimal length. Therefore we per- 
formed calculations with the computationally more ex- 
pensive algorithm, that allows for a change of the GS spin 
configuration. We considered system sizes L < 32 with 
up to 3000 samples and subsequent fits were restricted 
to 12 < L < 32. Albeit affected by finite size effects, we 
found d f = 1.080(5) (Q = 0.40) for the ±J-min2 DWs, 
compared to a value of df = 1.105(2) (Q = 0.43) when 
fitting the results for ±J-min MEDWs in the same in- 
terval. Hence, minimal-length and true minimum-length 
MEDWs are very similar, which means that the ± J-min 
MEDWs, where we can obtain results for large sizes, yield 
a reliable estimate of the behavior of shortest MEDWs. 

Also, we analyzed the scaling behavior of MEDWs of 
different energies for the ±J-modcl, in particular for 
those which have energies £qw = or Ebw = 2 (which 
constitute 96% of all MEDWs, the remaining 4% are 
Ebw = 4,6,8,10). In all cases we find again power- 



law behavior (not shown) with 



L f , the resulting 



exponents are almost compatible within error bars to the 
average result above. Note that in recent studies [13, |28| 
of the DW entropy, contrarily, the behavior of i?DW = 
and Ebw ^ DWs was very different. 

Also for the scaling of the domain- wall roughness with 
system size, we find a power law behavior, see Fig. [JJ] 
The resulting roughness exponents are shown again in 
Tab. IIIIl The roughness exponents, obtained using lo- 
cal slopes exhibit again larger error bars. Hence, in all 
cases studied here the roughness exponent appears to be 
compatible with unity, showing that the width of the 
MEDWs scales like the extension in y-dircction, as pre- 
dicted in Ref. [2lJ. Regarding the roughness, one finds 
stronger finite-size effects for small L, where we again 
have compared minimal-length and true minimum-length 
MEDWs. Fits to the data in the interval 16 < L < 32 
yields d r = 1.101(15) (Q ~ 0.75) for the true minimum- 
length ± J-min2 MEDWs and d r = 1.070(2) (Q = 0.79) for 
the minimal-length ± J-min MEDWs. Hence, minimal- 
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FIG. 6: (color online) Scale dependence of the average MEDW 
roughness. Fits to the form (r) ~ L dr for L > 50 (dashed 
lines) yield exponents listed in Tab. IIHI The inset shows the 
scaling of ±J-min compared to the ±J-min2 MEDWs. 



length and true minimum-length MEDWs are again very 
similar, which means that the ±J-min MEDWs yield a 
reliable estimate of the shortest MEDW behavior. 



is peaked close to (£). This holds also for MEDWs subject 
to ±J disorder (not shown), where one can observe an 
additional even/odd deviation in the distribution of the 
MEDW lengths. This deviation results in a preferential 
appearance of domain walls with even length, responsible 
for defect-energies £bw = mod(4). 

As it turns out, the particular shape of the scaling func- 
tion is very simple for the case of Gaussian disorder. It 
is possible to fit the distributions in this case satisfactory 
by use of log-normal scaling functions. 

Fig. [JJ shows an example, where the log-normal dis- 
tribution p(x) = exp [— (\\i(x/ n)) 2 /2a 2 ]/(xayf{2ii)) with 
x = £L~ d f was fit to the data (L = 160), resulting in a 
mean and standard deviation \x = 1.06(1) and er=0.24(l), 
respectively. This scaling function does not suit the bi- 
modal disorder case. Here, for MEDWs with minimal 
length, we observe an exponential decay of Pl{£) with 
increasing £ > (£) and estimated the decay exponent 
a = —5.9(2) from a fit to the data corresponding to 
L = 320. The DWs with maximal length somehow resem- 
ble the Gaussian case, but we could not find a meaningful 
distribution to describe it. 
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FIG. 7: (color online) Scaling plot of the MEDW length dis- 
tribution Pl{£) for Gaussian disorder, where the dashed line 
is a log-normal distribution with mean [i — 1.06(1) and stan- 
dard deviation a = 0.24(1), obtained from a fit to the data of 
L = 160. 



FIG. 8: (color online) Distribution Pl(w) of weights w on the 
MEDW segments for different system sizes L. The dashed 
line is a Gaussian distribution with mean n = —0.0003(40) 
and standard deviation a — 0.552(3), obtained from a fit to 
the data corresponding to L = 320. Results where averaged 
over 10 3 samples. 



So far, we have analyzed the scaling behavior of mean 
values, now we turn to the full distributions. The main 
result is that the distributions P L (£) of the MEDW length 
for different system sizes can be related to each other 
via a simple scaling relation. As can be seen from 
Fig. [7J rcscaling of the length distributions according 
to Pl{£) = L~ df f(£L~ df ) yields a collapse to a master 
curve. A qualitative similar behavior was found for the 
scaling of optimal paths in strong and weak disorder 
j39l . [4(j], the mass distribution of the backbone in crit- 
ical percolation [4l| and regarding undirected minimum- 
weight paths in 2d lattice graphs, where the effect of 
isotropically correlated bond weights on the scaling ex- 
ponents was investigated [42j ■ Note that the distribution 



Finally, we look at the distribution (w) of weights w 
which comprise the segments along the obtained MEDW 
for different system sizes, see Fig. [8] Since the segments 
available to the DWs are a result of a GS calculation, and 
because each MEDW is the result of another optimization 
procedure, the behavior is a priori not clear. Clearly, 
bonds with very negative energy will not occur, because 
these bonds with a large absolute value will be satisfied 
in the ground state with high probability, hence yield a 
positive contribution to the energy of a domain wall. This 
contributes to a decrease of the width of the distribution 
Pl{w), compared to the underlying disorder distribution. 
Apart from this effect, the resulting data compares well 
to a Gaussian with mean and a standard deviation close 
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to 0.5: from a fit to the data at L = 320, we have obtained 
a mean /j, = —0.0003(40) and a standard deviation a = 
0.552(3). Note that the distributions do not seem to 
depend strongly on the system size. 

IV. SUMMARY 

In summary, we have performed GS and DW calcu- 
lations for 2d Ising spin glasses with short ranged in- 
teractions via exact optimization algorithm. Exploring 
large system sizes, we investigated the fractal properties 
of MEDWs arising from Gaussian and bimodal disorder. 

Our approach is based on a minimum-weight path ap- 
proach for paths on undirected networks which can have 
negative edge weights. This allows for a more direct cal- 
culation of MEDWs. Presently, in the case of the de- 
generate ± J model, we are able to calculate the shortest 
and very long paths, which allows to obtain bounds on 
the fractal properties of the MEDWs. We believe that 
this approach could severe as a basis for the desired cal- 
culation of typical DWs for the ±J model. 

For the ±J model we found a lower bound df = 
1.095(2), which is clearly larger than unity, and we es- 
timated an upper bound df = 1.395(3) < 2. These expo- 
nents do not change when we analyze the scaling behavior 
of MEDWs restricted to certain energies. Our results for 
the length scaling exponent in case of Gaussian disor- 
der df = 1.274(2) is in agreement with prior results but 
with strongly enhanced precision. Further, it compares 
well to the result obtained from the recently proposed 
SLE scaling relation. Furthermore, we found that the 
full distributions of the MEDW lengths scale with the 
same fractal exponent df. This behavior was also found 



for observablcs in a different physical context. Finally, 
the width of the domain wall scales like its height, for all 
cases considered here. 

For later work, one could consider models with a T = 
transition between a ferromagnetic phase and a spin-glass 
phase, and study the behavior of the fractal dimensions 
around this transition. Here, from Eq. JT]), one would ex- 
pect to observe a discontinuous change of the exponents 
when approaching the transition, because 9 is expected 
to be the same everywhere in the spin-glass phase. Fur- 
thermore, it would be desirable to be able to sample GSs 
for the ±J model in equilibrium, hence each GS would 
contribute provably with equal weight/probability to all 
results. This would allow to determine a fractal dimen- 
sion of the average domain wall. This value might be 
related to the energy-scaling exponent 6 via a relation 
similar to Eq. (JTJ) . 
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